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ABSTRACT 

Recent X-ray and optical observations of the Perseus cluster indicate that a combination of 
weak shocks at small radii (> 20 kpc) and viscous and conductive dissipation of sound waves 
at larger radii is responsible for heating the intracluster medium and can balance radiative 
cooling of cluster cores. We discuss this mechanism more generally and show how the spe- 
cific heating and cooling rates vary with temperature and radius. It appears that this heat- 
ing mechanism is most effective above lO*^ K, which allows for radiative cooling to proceed 
within normal galaxy formation but stifles the growth of very massive galaxies. The scaling 
of the wavelength of sound waves with cluster temperature and feedback in the system are 
investigated. 
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1 INTRODUCTION 

The HQ-emitting filaments surrounding NGC 1275 in the core of 
the Perseus cluster (Lynds 1970; Conselice et al 2001) appear to act 
as streamlines revealing flows in the intracluster medium (Fabian et 
al 2003b). Many of the outer filaments are reasonably straight and 
stretch over 10s kpc. One dubbed the horseshoe (Conselice et al 
2001) folds back on itself just behind a large spherical-cap shaped 
depression in the X-ray intensity of the gas, suggestive of the flow 
on Earth seen behind rising spherical-cap bubbles in water (Fabian 
et al 2003b). Together these observations indicate that the intraclus- 
ter medium is not highly turbulent and is probably viscous, with a 
viscosity approaching the Spitzer-Braginskii (1956) value for an 
ionized plasma. 

Further evidence for the medium being viscous is provided 
by the shape of detached buoyant bubbles, seen most clearly in 
the Perseus cluster (Fabian et al 2003a). Simulations of bubbles 
rising in an intracluster medium show that viscosity enables bub- 
bles to remain intact for longer than a crossing time (Reynolds et 
al 2005). Strong magnetic fields are otherwise required to prevent 
bubbles from breaking up. 

If the intracluster medium is viscous then sound waves pro- 
duced for example by radio bubbles from the central active galaxy, 
as seen in the Perseus cluster (Fabian et al 2003a), can be rapidly 
dissipated. This provides an efficient mechanism for transferring 
energy produced by a massive black hole in a central cluster galaxy 
into the surrounding medium in a reasonably isotropic manner (see 
also Ruszkowski, Briiggen & Begelman 2004a,b). Although the en- 
ergy flow passes through highly directional jets, the resultant bub- 



bling produces approximately spherical sound waves which propa- 
gate into the cluster. 



Energetic outflows from black holes do not otherwise couple 
well to surrounding gas. Electromagnetic ones pass through the gas 
which, being ionized, is mostly transparent. CoUimated outflows, 
which tend to be the most energetic, are often too violent and bore 
right through the gas along the collimation direction (e.g. Cygnus 
A; Hydra A, McNamara et al 2005). Uncollimated winds might 
work but would tend to give strong shocks, which are not observed 
in the inner intracluster gas. Viscous dissipation can provide con- 
tinuous, gentle, distributed heat (Fabian et al 2003a), as is required 
(Voigt &. Fabian 2004). 



Here, we examine the dissipation of sound wave energy in in- 
tracluster and intragroup gas due to both viscosity and thermal con- 
duction. Both of these transport processes have a strong tempera- 
ture dependence which means that such heating dominates over ra- 
diative cooling above about 10^ K. Radiative cooling dominates 
below that temperature. This provides a clue as to why heating 
appears to stifle the cooling in clusters and groups (the cooling 
flow problem, see Fabian et al 2001) yet most galaxy formation, 
where much of the stellar component is a consequence of radia- 
tive cooling of baryons which have fallen into dark matter wells, 
has taken place. Indeed the process which stifles cooling in mas- 
sive systems probably is responsible for truncating the upper mass 
range of galaxies (Fabian et al 2002; Benson et al 2003; Binney 
2004). 
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2 VISCOUS AND CONDUCTIVE HEATING 

2.1 Sound wave heating and radiative cooling rates 

We now consider the ICM heating rate caused by the dissipation of 
radio-galaxy induced sound waves. Suppose the ICM possesses a 
kinematic viscosity i/ and a thermal conductivity of k. The absorp- 
tion coefficient for the passage of sounds waves is given by (Landau 
& Lifshitz 1987) 



7 = 



4 ft / 1 
— 

3 p \cv 



(1) 



where / is the frequency of the sound waves, Cg is the sound speed, 
p is the density of the ICM, and cv and Cp are the specific heats at 
constant volume and pressure, respectively. The dissipation length 
(i.e., the length scale over which the energy flux of the sound wave 
decrease by a factor of e) is given by £ = I/27. Noting that 
Cp/cv = 5/3 for a fully ionized gas, we can write 



£ = 



3c? 



Stv^P [2v + K/pcp] ■ 



(2) 



We note that this expression is only strictly true in the regime in 
which the wavelength of the sound waves are short compared with 
the dissipation length I. We briefly discuss a possible violation of 
this condition later in this Section. Now, suppose that the viscosity 
and thermal conductivity are fixed fractions, ^„ and of their 
unmagnetized values. For a fully ionized hydrogen plasma, in c.g.s. 
units, we get 



1.0 X lO'^^T^^^n' 
2.36 X 



(3) 
(4) 



pCp 

where T7 = T/(10'^K), and n is the electron number density, 
and we have taken the Coulomb logarithm to have a value of 37. 
Substituting the numerical values into eqn.|2| we get 



: 697- 



nTi\f: 



0.1 J 



+ 11.8 



k 0.1 , 



■ kpc 



(5) 



where we have expressed the sound wave frequency in units of 
per megayear (i.e., f-e = f / iW~^ yr~^) and have normalized 
to plausible values of viscosity and thermal conductivity (follow- 
ing the arguments of Narayan & Medvedev 2001). Hereafter, we 
shall denote 



^ 0.1 0.1 



(6) 



We note that, since temperature and density vary with radius in the 
cluster, the dissipative effects of viscosity and thermal conductiv- 
ity vary but maintain a fixed ratio (as a consequence of assuming 
that and are constant). We also note that thermal conduction 
is over an order of magnitude more effective at dissipating sound 
wave energy if ~ ^k. 

We now model the central AGN as a source of acoustic energy 
at the center of the cluster. Suppose that the AGN injects an acoustic 



luminosity of Linj into the ICM at an inner radius of r = 



If 



Ls{r) is the acoustic luminosity at radius r > rin, the definition of 
dissipation length gives 



dLs _ 
dr 

which has the solution 

Ls (r) = Linj exp 



■ dr 



(7) 
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Figure 1. Dissipative volume heating rates (solid lines) and radiative vol- 
ume cooling rate (dashed line), using the temperature and density profile of 
(Top) A 2199 (Johnstone et al. 2002) and (Bottom) the core of the Perseus 
cluster (Sanders et al 2004). The heavy solid line is for the following set 
of illustrative set of parameters for which the heating and cooling approxi- 
mately balance; L44 = 1 for A2199, 5 for Perseus, /—g = 0.1, £,1, = 0.1, 
^k; = 0. Other choices of the acoustic power (Z/44 will simply rescale 
this line in the vertical direction. Other choices of /_6, = 0.1, and 
= change the dissipation lengthscale as well — the family of ad- 
ditional (thin) solid lines shows the effect of successive doublings of fre- 
quency (/-6 = 0.2, /_6 = 0.4, and /_6 = 0.8). The vertical dashed 
lines indicate a fiducial cooling radius where the radiative coohng time of 
the gas is 3 Gyr 



The volume heating rate due to viscous and conductive dissipation 
is 



Ls{r) 



which, for r <^ lis approximately 

2TvfLinj [2U + K/pCp] 



Edis 



and can be evaluated to give 



3ch 



(9) 



(10) 



ediss « 3.8 X 10-2^T7n-V-6r2"'i44Cergcm-='s-i , (II) 

where L44 = Linj/ (10''* ergs"' ) and r2 = r /(ICQ kpc). 

As illustrations of the radial dependence of the viscous and 
conductive dissipation, we compute ediss using the temperature 
and density profiles of Abell 2199 from lohnstone et al. (2002; 
T = 4.4r§-^ keV and 71 = 6.0 x IQ-V^" '^^ cm"^ ) and beyond 
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about 20 kpc in the core of the Perseus cluster from Sanders et al 



(2004; r 



b.MrV kcV and n 



^cm-^). 



We compare these with the radiative cooling as derived from the 
formula of Tozzi & Norman (2001) which, expressed in our units, 
becomes 

£,ad = IQ-^^mne [l.lSTf + 5.3T7°-^ + 6.3] (12) 

As shown in Fig. 1, the dissipative heating can be comparable or 
exceed the radiative cooling within the fiducial cooling radius in 
which the cooling time r = 3nfcT/2ei.ad exceeds 3 Gyr. Indeed, in 
the case of Abell 2199, the dissipative heating closely balances the 
radiative cooling out to the cooling radius for the following choice 
of parameters; L44 = 1, / = 0.2, = 0.1, ^„ = 0. The dis- 
sipation lengthscale is approximately 100 kpc, and the heating rate 
displays power-law radial profile within this radius (reflecting the 
simple power-law forms of the assumed temperature and density). 
A sharp (exponential) cutoff in the heating rate is seen for radii 
larger than the dissipation length. Varying the acoustic luminosity 
of the AGN simply changes the normalization of this heating rate. 
On the other hand, varying the frequency of the sound waves leads 
to dramatic changes in the dissipation lengthscale. In particular, in- 
creasing the sound frequency while leaving the other parameters 
at the values listed above rapidly brings the dissipation lengthscale 
within the cooling radius. A detailed balance of heating and cooling 
is then no longer possible — the innermost regions can be strongly 
heated while radiative cooling occurs unchecked at the cooling ra- 
dius. In the case of the Perseus cluster, the steep density profile re- 
sults in the radiative cooling possessing a significantly stronger ra- 
dial dependence than the "monochromatic" (i.e., single frequency) 
dissipative heating profile. Thus, any balance between sound wave 
dissipation and radiative cooling must employ acoustic waves with 
a range of frequencies. 

To quantify this, suppose that the injected acoustic luminosity 
between frequencies / and / + d/ is Linj (/) d/. Furthermore, sup- 
pose that all frequencies are injected at the same radius r = nn. 
The full expression for the dissipative heating is then. 



edii 



Jo 



exp 



dr df. 



(13) 



For comparison with Perseus, it is useful to quantify the flattening 
of the dissipation profile due to the spectrum of frequencies in the 
case of power-law radial dependences of temperature and density. 
Suppose that T oc r" and n oc r~'^, where we expect a, /? > on 
general physical grounds. Using eqn (11) we can then approximate 
the dissipative heating profile of a given frequency as a power-law 
in radius. 



ediss(/) oc r 



(14) 



which is truncated at radius r = I. Suppose that the spectrum 
of frequencies has a power-law form, Linj{f) oc /^"' between 
/ ~ /min and / = /max. Noting that £ increases with decreas- 
ing frequency, we can delineate the following regimes; 

(i) If r > ^(/min), the radius under consideration is outside of 
the overall dissipation region and the dissipative heating profile is 
rapidly declining. 

(ii) If r < i'(/max), even the highest frequency component has 
not yet dissipated away and the radial dependence of the heating 
is just the sum of similar monochromatic components, ediss oc 

(iii) If £(/max) < r < ^(/min). Only the frequency components 
upto f = f contribute, where / is defined by r = £(/). In this 



Cdiss OC 



^f-'df 



(15) 



For 7 > 3, the lowest frequencies always dominate the energetics 



and we again get the monochromatic case eaiss oc 
7 < 3, we get £diss oc , where 



_c«+/3-2 



For 



(16) 



The second term on the RHS of the last expression represents the 
modification of the monochromatic result due to the spectrum of 
frequencies. 

Using these results, we can predict that for dissipative heating to 
balance radiative cooling in the Perseus cluster, we require a spec- 
trum of fluctuations with 7 « 1.8 across at least the frequency 
range /_6 = 0.2 - 1. 

For a given suppression factor (^^ and ^k), thermal conduc- 
tion can be significantly more important than viscosity at dissipat- 
ing sound wave energy. However, this is only strictly true when the 
dissipation wavelength is large compared with the wavelength of 
sound. If the thermal conduction is such that the dissipation length- 
scale I is comparable to or smaller than the wavelength of the sound 
wave under consideration, dissipation due to thermal conduction 
will be ineffective and the sounds wave will propagate as an isother- 
mal sound wave. Indeed, if thermal conduction occurs at the level 
suggested by Narayan & Medvedev (2001), 

~ 0.2, (17) 

the dissipation lengthscale would be extremely short for all but 



the lowest frequency soundwaves 



150/_g pc using n 



0.02 cm~ and T7 — 4). Thus an implicit prediction of Narayan & 
Medvedev (2001) is that the sound waves observed in the Perseus 
cluster are isothermal sound waves. Indeed, the thermal conduc- 
tion can become so ineffective at dissipating acoustic energy that 
viscous dissipation dominates even if > C^- ■ 

The theory outlined above only strictly applies to the case 
of linear sound waves; this imposes a radius dependent maximum 
acoustic luminosity to which this theory can be applied. Expressed 
in terms of pressure fluctuations, 6p, the acoustic luminosity cross- 
ing a shell at radius r is, 

i5p' 



4-Kr 



(18) 



The acoustic waves become strongly non-linear when Sp w p. 
Thus, the maximum acoustic luminosity for which our model ap- 
plies is 

2 

Ls,max « 47rr^^ fa 3.1 x 10"*^ nTy' Vi ergs'^ . (19) 

Using this equation, we see that the fiducial model described above 
for A 2199 (L44 = 1, / = 0.1, = 0.1, = 0) is only valid for 
radii greater than about 20 kpc. Within that radius, the energy must 
be transported within the form of shocks and/or strong convective 
motions. 

For comparison, we now consider the dissipation rate in a 
weak shock, tshock oc {6P/P)^pf (Landau & Lifshitz 1987; Stein 
& Schwartz 1972) 



Cshock ' 



4 X 10"^* Lffn-^/^T, 



1/2^-5/4 



/-6 erg cm 



(20) 

This is comparable with the required heating rate at about 10- 
20 kpc radius in the Perseus cluster (where a weak shock is seen. 
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Fabian et al 2003a) if /g ~ 0.1 and L44 ~ 10. For the den- 
sity and temperature profiles of that cluster, ediss cc r^'^'^" and 
Cshock oc r^'^. Viscous dissipation is therefore most effective at 
larger radii and for higher frequencies, shock dissipation applying 
most to small radii and low frequencies. Since 96 per cent of the 
volume within 60 kpc lies beyond 20 kpc, it is reasonable to state 
that viscous heating can dominate the volume heating. 

2.2 A note on the effect of magnetic fields in the ICM 

In the above analysis, the viscosity and thermal conductivity were 
assumed to be a fixed fraction of the Spitzer-Braginskii (1956) 
value for an ionized gas. Cluster cores do have significant mag- 
netic fields (e.g. Carilli & Taylor 2002) which may of course alter 
the viscosity. What effect this has is unclear. 

There is, however, evidence from Faraday rotation measures 
(RMs) observations (Carilli & Taylor 2002 and references therein) 
that the field is coherent in patches larger than the ion mean free 
path. Typical values of the cell size estimated from RM distribu- 
tions and quoted in the literature are 5-10 kpc. A power-spectrum 
analysis of Hydra A, a bright radio galaxy and one of the best cases 
for probing a cooling core cluster as it provides a probe of the RM 
distribution on scales from 0.2 to 40 kpc, has been carried out by 
Vogt & Esslin (2005), who find a magnetic autocorrelation length 
of 3 ± 0.5 kpc. This autocorrelation length can be compared to 
the ion mean free path (mfp) of Q.2ZT^n~^ kpc for typical clus- 
ter core parameters of Tc = 3 x 10^ K and n^T^ = 0.01 cm^'^ . 
Since the mfp is smaller than the dominant scale size of the mag- 
netic field, the viscosity along the magnetic field direction will be 
high for much of the volimie under consideration. This is also rele- 
vant to the issue of conduction, which can also dissipate the energy 
in sound waves (Landau & Lifshitz 1987; sec also Ruszkowski, 
Briiggen & Hallman 2005), but is not strong enough to offset cool- 
ing in many clusters by itself (Voigt & Fabian 2004; Kaastra et al 
2004). 

The simple model considered in Section 2. 1 can be motivated 
from the notion that intracluster medium consists of many magnet- 
ically isolated blobs of gas, each of size ranging over a few kpc. 
Within a blob the field is fairly coherent. Conduction of heat from 
the outside of the core to the centre is inhibited by the boundaries 
between the blobs. The energy in sound waves propagates freely 
across the blob boundaries and can be dissipated within the blobs 
by both conduction and viscosity. 

There are other possible sources of viscosity in a tangled mag- 
netic field. For example if the field is in the form of flux tubes there 
can be drag between the high field regions and the low field ones 
(Longcope, McLeish & Fisher 2003). Alternatively, the spatial con- 
figuration of the field may change as sound waves pass across a 
region in the sense that shaking up the field causes it to shift to a 
lower energy configuration, so releasing energy either immediately 
or later. This would operate as a bulk viscosity. (The issue of bulk 
viscosity is also mentioned by Ruszkowski et al 2004). 



3 THE WAVELENGTH AND FREQUENCY OF SOUND 
WAVES 

Subcluster mergers are conmion in clusters and these will generate 
large pressure discontinuities and sound waves. The wavelengths 
of sound waves remaining after shocks have passed will presum- 
ably be on the scale of cluster cores or hundreds of kpc. On these 
scales the gas may be turbulent so the motions may degrade down 



to scales of tens of kpc where viscous dissipation can operate. The 
outer gas in a cluster may therefore be noisy with a wide spectrum 
of wavelengths. Pringle (1989) noted that the dense regions of a 
cluster core can focus sound waves inward, meaning that signifi- 
cant heating may result at the centre. Fujita et al (2004), using a 2D 
simulation, claim that such motions generate strong turbulence. Es- 
tablishing the energetics and efficiency of such processes requires 
numerical calculations beyond the scope of this paper. We note that 
propagating turbulence inward, where the gravitational potential is 
steep and the gas dense, is energetically challenging. 

A central active nucleus giving rise to jets will inflate cav- 
ities of relativistic radio-emitting plasma within the intracluster 
medium. Such cavities or bubbles are seen in the Perseus cluster 
(Bohringer et al 1993; Fabian et al 2000; Fabian et al 2003a), Virgo 
cluster (Wilson et al 2002; Forman et al 2003), A2052 (Blanton 
et al 2001, 2003), A2597 (Nulsen et al 2002), A4059 (Heinz et al 
2002) and Hydra A (McNamara et al 2000). Such bubbles inflate 
and then break away due to buoyancy (Churazov et al 2000, 2001) 
and so even a constant jet will lead to a bubbling effect. This gen- 
erates a repetitive pressure excess and thus a sotmd wave with a 
period equal to the repetition time of the bubbles. 

A bubble expands because of its excess pressure relative to 
the intracluster medium at a rate which varies as t^^^, where t is 
time. It detaches from the jet when buoyancy dominates which is 
roughly when its expansion velocity drops to about one half of 
the local Keplerian velocity due to the local gravitational field of 
the central galaxy and cluster, vk (Churazov et al 2000). A new 
bubble then begins to grow. This gives a natural bubbling period 
', where P is the pressure of the intracluster 
gas and L is the power of the jet. Now the central cooling time of 
many clusters with cool cores levels off at about 10* yr within a ra- 
dius of 10 kpc (Voigt & Fabian 2004) which implies that (between 
clusters) P cx T'""^^ . If wc also assume that vk tx T^'^ where the 
gas temperature is T. Then / oc ^^-1/2^3/2 fvuther scaling 
depends on L. 

The dependence of L on T is probably fairly steep. If we 
assume that the jet power L does balance cooling then it is pro- 
portional to the X-ray luminosity Lcooi, within the cooling radius 
(where the cooling time equals the cluster age). Lcooi is tabulated 
for nearby clusters by Peres et al (1998) and using these values and 
the cluster temperatures we find (Fig. 2) that L oc T" with a ~ 4. 
Therefore / oc T~^^^ . As a rough check, the scaling predicts that 
radius R when a bubble separates, R oc L''^P-^I\-"^ oc T. 
This scaling result is shown in Fig. 3 with results for individual 
clusters drawn from the compilation of Dunn & Fabian (2004). 

The luminosity drops out of the dissipation formula (11) 
which now varies steeply with temperature; eaiss oc T'^/'^r~^ and 
at a given radius ediss/frad oc T^. 

Although we have in this Section identified the bubbling fre- 
quency as the dominant one, we speculate that the underlying vari- 
ability of the radio jets occurs at higher frequencies and carries 
much power. Most jets are not constant as shown by repeated out- 
bursts on small scales (e.g., in Perseus A=3C84, Taylor & Ver- 
meulen 1996, see also Reynolds & Begelman 1997, and in general 
Zensus 1997). These variations will be transmitted through to the 
surface of the bubble leading to its growth being erratic and a spec- 
trum of high frequency sound waves being present at small radii, 
with the bubbling frequency being the low-frcqucncy cutoff. The 
steep dependence on frequency means that higher frequency waves 
are damped most easily. A simple estimate of the heating effect of 
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Figure 2. The X-ray luminosity witliin the cooling radius plotted against 
the cluster temperature, using results for cooling flow clusters taken from 
Peres et al (1998). The best fitting power law is shown, with index 3.8. 
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Figure 4. Heating (sohd) and cooling (dotted) rates plotted as a function of 
temperature for a pressure nT = 10® cm~^ K and radius of 10 kpc. For 
the heating rate, f^^Lnl = 1. 
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Figure 3. The radius of bubbles is shown as a function of the surrounding 
temperature. Data are for attached bubbles from the compilation of Dimn & 
Fabian (2004). The open point is for A2052 where there appears to be much 
cool gas accumulated around the rims (Blanton et al 2003) so the points 
may reasonably shift to the right. The dashed line is not a fit but indicates 
the /t oc T relation tentatively derived in the text for the rough radius, R, 
at which bubbles detach from the nucleus. Bubbles should therefore only 
exist below, or close to, the line. 



bubbles from the bubble size and recurrence rate (e.g. Birzan et al 
2004) will underestimate the true heating rate. 

As an indication of the temperature dependence of the dissipa- 
tion versus radialivc cooling we show Cdiss and trad as functions of 
temperature in Fig. 4. We assume that f-eL44^ = l,r = 10 kpc 
and P = nT = 10^ cm~^ K which means that if the above scal- 
ing is correct, the heating rate in real objects is steeper. The plot 
shows that sound wave dissipative heating dominates above about 
10^ K and radiation dominates below. The process is therefore rel- 



evant as an explanation for the truncation of cooling, and thus of 
star formation and the total stellar mass, in massive galaxies. 



4 DISCUSSION 

We have shown that viscous heating of intracluster or intragroup 
gas can effectively compete with radiative cooling above gas tem- 
peratures of about lO'^ K. The repetitive PdV work done by the 
formation of bubbles is dissipated by weak inner shocks and outer 
viscous processes. More rapid intrinsic variability in the jets leads 
to quasi-continuous generation of higher frequency sound waves 
which dissipate throughout the cooling region. Sound waves pro- 
duced by jets from a central black hole can thereby strongly reduce 
cooUng in the surrounding gas and prevent further growth of the 
central galaxy by accretion of radiatively cooled gas. This mecha- 
nism may explain the upper-mass cutoff in galaxies since it operates 
best at high temperatures. 

If cooling initially dominated in a cluster core then the tem- 
perature profile would have dropped to the local vrrial one, which 
for a NFW potential (Navarro, Frenk & White 1997) approaches 
T oc r''^'^ for the central cusp (i.e. within the central 100 kpc). 
Steady heating then applied to this gas will tend to push the gas out- 
ward, lowering its density whereas continued cooling means that 
gas flows inward, in both circumstances following close to the lo- 
cal virial temperature profile. This roughly explains the observed 
temperature profiles, given that much of the inner few kpc are oc- 
cupied by bubbles. 

The prevalence of the central X-ray surface brightness peak 
characteristic of short radiative cooling times (Peres et al 1998; 
Bauer et al 2005) and the remarkable similarity of the cooling time 
profiles in many clusters (Voigt & Fabian 2004) suggests that a 
good heating / cooling balance is commonly established. This pre- 
sumably requires feedback with the central power source, the ac- 
creting black hole. How that is established so that the properties of 
the gas at say 50 kpc are attuned to those at the black hole accretion 
radius of say 100 pc, is not completely clear. 

A partial answer can be seen from the following. If the radio 
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source becomes too vigorous, then most of its power will be de- 
posited at large radii beyond the cooling radius both if it becomes 
an FRII (as in Cygnus A) or remains an FRI when the associated gi- 
ant bubbles (as in Hydra A, Nulsen et al 2004, or MS0735. 6+7421, 
McNamara et al 2004) create low frequency global disturbances 
which dissipate poorly. If the radio source shuts off due to lack of 
fuel, then the inner gas will cool and, angular momentum permit- 
ting, accrete and provide that fuel. Rapid variations in jet power will 
lead to high frequency disturbances which dissipate and so rapidly 
heat those inner regions. A trickle of cooling and accretion explains 
the prevalence of warm optical nebulosity (Crawford et al 1999 and 
references therein) and cold gas (Edge 1991) near the centres of 
these objects. The inner cooling time of these objects at ~ 10* yr 
(Voigt & Fabian 2004) is similar to the dynamical timescale of the 
central galaxy. 

Continued heating by a central AGN in cluster and group cores 
means that the central black hole continues to grow. The total ac- 
creted mass can be considerable (see Fabian et al 2002; Fujita & 
Reiprich 2004). This may skew the mass - velocity dispersion re- 
lation (M — a) for massive black holes upward at the highest black 
hole masses, such as expected of the massive black holes in groups 
and clusters. 

Finally, we note that a viscous intracluster medium means that 
motions in the cores of relaxed clusters will be damped and, apart 
from the regions directly around the inflating and rising bubbles, or 
where the dark matter potential or galaxy are sloshing around, fluid 
motions should be small and highly subsonic (< 0.3cs). Of rele- 
vance here is the overall scale of any turbulent energy cascade, i.e. 
the difference in length scale between the energy injection scale, 
presumably the bubble size, and the dissipation scale, which may 
be as large as the size of magnetic cells or similar distinct struc- 
tures or extend down to much smaller scales. Measurements of the 
widths and shape of iron-K emission-lines with Suzaku will help 
to distinguish laminar motions from turbulence (see Inogamov & 
Sunyaev 2003 for predictions of turbulent profiles). 
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